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This manuscript reports research performed for NASA-MSFC during the period 
from February 1 to December 31, 1991 under Contract #NAS8-36955-114. The objective of 
this project was to investigate the feasibility of using frequency response techniques for 
enhancing destructive physical analysis and for nondestructive testing of aerospace bat- 
tery electrodes. Nickel, cadmium, silver and zinc electrodes were tested by imposing al- 
ternating current upon the electrodes and measuring the magnitude and phase of the 
response voltage. This yields an impedance spectrum for the battery electrode from 
which electrochemical kinetic, double layer capacitance and mass transfer effects can be 
characterized. Frequencies from 10 kHz to 0.1 mHz were used in the testing. 


Experimental 

Figure 1 shows a diagram of the experimental apparatus. The centerpiece is a 
Schlumberger Model 1255 Frequency Response Analyser (FRA) and Model 1250 Electro- 
chemical Interface (El). The FRA contains a signal generator , the output of which can 
modulate the voltage or current output of the El. The FRA also has two signal inputs 
which are connected to the two analog outputs of the El, proportional to the cell current 
and the electrode potential. 

The battery electrodes were placed in prismatic cells fabricated from plexiglass. 

These cells typically contained one working electrode and one or two counterelectrodes. 

The working electrode potential was measured relative to a reference electrode that was 
also in the cell. Only a small amount of current, on the order of microamperes, is drawn 
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from the reference electrode during potential measurements, hence the working electrode 
polarization can be measured without interference from the counterelectrode. A silver 
reference electrode was used for silver and zinc working electrode measurements while a 
nickel reference electrode was used for nickel and cadmium electrode measurements. 

The FRA and El were controlled by an Apple II computer through a GPIB inter- 
face. Programs written in Applesoft Basic were used to control the amplitude and fre- 
quency during logarithmic sweeps of the working electrode impedance versus frequency, 
and to control the El during charge and discharge of the cells. The results were displayed 
on the computer screen and recorded in a notebook for processing. Bode magnitude plots 
(impedance magnitude vs. log frequency) and Bode angle plots (impedance phase angle vs. log 
frequency) were used to represent the data. 

Frequency response analysis is valid only for linear systems, but electrochemical parame- 
ters such as interfacial resistance, double-layer capacitance and diffusion (Warburg) impedance are 
strong functions of potential. Hence, it is important that the amplitude of the voltage response be 
small. Each of the above named impedances can be treated as constant if the voltage amplitude 
driving current through that impedance does not exceed 5 mV. The computer programs were writ- 
ten so that the total working electrode voltage amplitude was about 5 mV so that the amplitude 
through the individual impedances would be no larger than this maximum permittable value. 

Mathematical Modeling 

This section will briefly summarize the modeling techniques for analyzing impedance data. 
More complete presentations can be found in the literature. Two techniques are reported- equiva- 


MacDonald, 1991). Two techniques are reported- equivalent circuit analysis and finite difference 
numerical analysis. Most of the work in this project used the first technique, but work with the 
second technique was begun. 

Equivalent Circuit Analysis- The traditional approach to electrochemical impedance data is 
to model it in terms of equivalent electrical circuit elements, fitting the experimental data to the cir- 
cuit with the method of least squares. Equivalent circuit elements used in this work will be dis- 
cussed briefly. 

The double layer capacitance (C ) of the electrode/solution interface is usually modeled as 
a pure capacitance. The double layer itself is a solution layer.about 10 to 20 Angstroms thick, adja- 
cent to the electrode surface which can be charged due to: 1) nonrandom alignment of dipoles 
near or at the electrode surface or 2) preferential adsorption of either the anion or cation at the elecr 
trade surface. This capacitance is relatively large, on the order of 10' 5 F/cm 2 , and is strongly poten- 
tial dependent, but can be treated as constant for potential variations of the order of 5 mV or less 
(Newman, 1991). 

There are two electrochemical phenomena that can be modeled as pure resistance with no 
reactive component: ohmic electrolyte resistance and Faradaic resistance. The ohmic electrolyte 
resistance R q is usually linear over the entire range of potential and is simple to analyze. In the 
case of alkaline aerospace batteries the electrolyte conductivity is quite high so that R q is low. The 
Faradaic resistance R is the resistance to charge transfer across the electrode/electrolyte interface 

P 

and is also known as the kinetic resistance. It is highly non-linear since the relationship between 
electrode reaction rate (current density) and interfacial potential difference is exponential, but the 
exponential relationship can be linearized for potential variations of about 5 mV or less. 



The Warburg impedance is caused by concentration gradients in the mass transfer bound- 
ary layer (thickness of 100 |im order of magnitude) adjacent to an electrode, has both resistive and 
reactive components, and increases as frequency decreases. This impedance becomes appreciable 
when the frequency is low enough that significant depletion of reactants or accumulation of prod- 
ucts occurs during the anodic or cathodic half of a sinusoid. The Warburg impedance for an elec- 
trode in a large excess (semi-infinite) of solution is 

Z = a/(jcoD) 0J (1) 

w 

where D is the diffusion coefficient, co is the frequency and a is the Warburg coefficient, 
a=RT/C (nFf (2) 

b 

R is the gas constant, T the absolute temperature, C the bulk concentration, n the stoiciometric 
number of electrons transferred and F is Faraday's constant (96487 Coulombs/mole electrons). 

The phase of the Warburg impedance is -45° and the magnitude is proportional to the inverse 
square root of frequency. More complex expressions are required when there is not excess solu- 
tion or when there is more than one diffusing species contributing to the impedance. An exponen- 
tial relation between concentration and potential (the Nemst equation or something similar) means 
that this impedance is also highly non-linear but can be linearized for voltage amplitudes less than 
5 mV, 


Another circuit element that has gained popularity in recent years,especially when the data 
to be analyzed does not fit simple circuit elements, is the constant phase element (CPE). It has the 
mathematical form 


Z =A(jw)“ 


cpe 


( 3 ) 


where A and n are constants that can be fit to experimental data. The impedances R q , R , C and 



Z can be thought of as CPE's where m is a theoretical valiie. The CPE can be used to account for 

W 

complex phenomena such as coupled diffusion effects or porosity fluctuations. It is often used as 
an empirical fitting device and has the disadvantages associated with such tools. However, if the 
alternative is an equivalent circuit with an excessive number of elements (circuits with 15 or 20 el- 
ements have been proposed in some systems without physical explanation of the elements) then 
the CPE may be an attractive alternative. 

Figure 2a shows a simple equivalent circuit that has been used often. C is considered par- 
allel to R and Z , forming an equivalent interfacial impedance that is in series with the R .. The 

f w Q 

justification cited is that the first three can only be determined by measurements at an interface. 
However, since R and C are double-layer phenomena while Z is a diffusion layer phenomena 
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(orders of magnitude thicker than the double-layer), the circuit of Figure 2b seems more physically- 
realistic. 

The two preceding equivalent circuits are for smooth planar electrodes. Battery electrodes 
usually are porous electrodes in which there is pore electrolyte resistance as well as the interfacial 
impedances discussed in the preceding paragraphs. Figure 2c shows transmission line model that 
is often used to model porous electrodes, where R q is the pore electrolyte resistance and Z repre- 
sents the collective interfacial impedances. The impedance of this electrode, when there is no D.C. 
bias current and no concentration gradients along the pore axis (DeLevie, 1963), is 
Z = cosh(mL)/(K msinh(mL)) (4) 

p eff 

where 

m 2 = a/K Z, (5) 

eff 1 

a is the interfacial surface area per unit volume and k is the effective conductivity of the elec- 
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Figure 2. Equivalent Circuits used in Impedance 
Modeling. 






trolyte. In the limit of small mL 

Z = Z/aL (6) 

P 1 

while in the limit of large mL 

Z = 1/K m (7) 

p eff 

At this latter limit, the phase angle of Z is halved so that double layer capacitance has an apparent 
phase angle of -45° while semi-infinite Warburg impedance has an angle of -22.5°. 

The equivalent circuit for Z used most in this study is shown in Figure 2d, and is similar to 
Figure 2b except that the Warburg impedance was replaced by a CPEand the ohmic resistance 
was deleted. Also, almost all the phase angles measured were less than zero (capacitive reac- 
tance), and when it is said that a phase angle increases it is referring to the magnitude of the angle. 

Numerical Finite Difference Models- The physics and chemistry that govern battery elec- 
trode operation can be expressed in the form of coupled differential equations with much greater 
clarity and flexibility than is possible using equivalent circuit techniques. In particular, compli- 
cated phenomena that would otherwise get lumped into CPE's can be modeled in this way. Nu- 
merical finite difference solutions for models of charging or discharging batteries abound in the 
literature. However, these techniques have not been widely utilized to analyze impedance data. 
Results and Discussion 

Silver/Zinc Cells- A cell consisting of one silver and one zinc electrode from a Yardney 50 
Ahr battery was assembled and tested. The reference electrode was an anodized silver wire that 
was connected to the silver working electrode during cycling. The cell was cycled twice, with dis- 
charge at 100% of measured capacity, and the impedance spectrum of each electrode was mea- 
sured at 100%, 75%, 50%, 25%, and 0% states of charge (SOC). 



Figures 3 and 4 respectively are the Bode magnitude and angle plots for the silver electrode 
during discharge on the second cycle. It shows a trend of decreasing magnitude and phase angle 
with increasing state of charge at low frequency, reflecting increased mass transfer impedance, al- 
though most of the lines cross at 3 mHz. These general trends were observed in both cycles. 
Figures 5 and 6 are the plots for the zinc electrode during discharge on the second cycle. In this 
case the fully charged state has the highest impedance magnitude and phase angle. In both elec- 
trodes it seems that the bare metal in contact with solution has the higher impedance, while the 
metal with the fully oxidized surface has the lowest impedance. 

Tests were also made with two Yardney 150 Ahr cells. These cells had been cycles at least 
100 times, and as received had rest potentials of 1.6V (cell #1) and 0.3 V(cell #2). Impedance 
measurements were first made for the two cells as received and shown in Figures 7 and 8. Cell 
number two had a much higher low-frequency impedance. It is thought that the zinc electrode 
morphology in this cell had deteriorated, reducing the zinc active surface area and increasing the 
mass transfer limitations. The cells were then charged to 1.8 V rest potential and spectra were 
again run and are shown in Figures 9 and 10. Charging increased the magnitude of the impedance 
for both cells but caused the phase angles to become less negative at low frequency. The work re- 
ported in the previous paragraph showed that impedance magnitude is highest for bare metal elec- 
trodes (the zinc electrode in this case), hence the zinc electrode impedance increases more than the 
silver electrode impedance decreases during charge. It is also interesting to note that the charged 
cell has positive phase angles at high frequency. It has been shown in other systems(Keddam, et 
al., 1981, 1984; MacDonald, 1991) that this can result from non-elementary electrochemical 
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reactions. This most likely occurs at the zinc electrode in which the reaction product, zincate ion, 
is quite soluble in the electrolyte solution. 

Figures 1 1 and 12 compare the experimental results for the impedance of a discharged 
silver electrode with the best fit to the equivalent shown in Figure 2d at the limit given by Equation 
6. The circuit seems capable of following the trends in the magnitude, while the maximum angle 
discepancy is about 0.25°. Table 1 gives values for the equivalent circuit parameters. The 
parameter aC, which is the capacitance per unit volume, is a figure of merit that characterizes the 
electrochemically active surface area. 

Nickel! Cadmium Cells- Impedance spectroscopy was used by Armstrong et al. (1975) and 
Armstrong and Edmondson (1974) to study cadmium electrodes while Lenhart et al. (1988) have 
studied porous nickel electrodes. NiCd cell studies have been performed by Sathyanarayana et al. • 
(1979) and Zimmerman et al. (1982). 

A cell consisting of a nickel electrode sandwiched between two cadmium electrodes was 
constructed, and impedance spectra for the nickel and cadmium electrodes were obtained. Figures 
13 and 14 compare the spectra for charged and discharged cadmium electrodes. As with the zinc 
and silver electrodes, the low frequency impedance magnitude and phase angles for the bare metal 
(charged cadmium) are higher than for the oxidized metal. There was little observed dependance 
of impedance with state-of -charge for the nickel electrode, which is always in an oxidized state. 

Figures 15 and 16 compare the cadmium electrode experimental impedance with the 
(Figure 2d) equivalent circuit best fit Of all the electrodes studied, this one seemed to fit the 
equivalent circuit best. Figures 17 and 18 compare the results for the nickel electrodes, which give 
a fairly reasonable agreement between experiment and equivalent circuit. There is some 
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discrepancy in the angle plot although the trends are comparable. The results are markedly 
different from Lenhart et al (1988), who reported all positive phase angles for porous nickel 
electrodes over approximately the same frequency range. This may be due to the presence of 
LiOH in the electrolyte and to differences in electrode manufacture. 

It should be emphasized that the fit of equivalent circuits to experiments were obtained by 
replacing the Warburg impedance, which is based on fundamental electrochemical transport 
principles, with the empirical constant phase element. Replacement of equivalent circuit models 
with coupled differential equation models is recommended for future studies. While this approach 
is much more difficult to implement it is expected to provide a more rigorous understanding of 
porous battery electrodes. 

It is also interesting to note the wide variation in a C values in Table 1. It is possible that 
the CPE is also accounting for the some or most of the double-layer capacitance effects, and that a 
more fundamental modeling approach will yield double-layer capacitances that do not differ by so 
many orders of magnitude. However, other workers have also observed order of magnitude 
differences in capacitance effects between electrodes. For instance, Tiedemann and Newman 
(1976), using models without empirical components, observed that capacitance effects for porous 
lead electrodes are two orders of magnitude higher than for porous lead dioxide electrodes. This is 
because double layer capacitance is a chemical, and not purely an electrical, effect. 

Conclusions 

1. State of charge estimations can be made using electrochemical impedance techniques. 

In particular, the low-frequency impedance of metal electrodes is higher in the fully reduced state 
than in the fully oxidized state. 
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2. Comparisons between healthy and unheajthy AgZn cells show that, at least in some 
instances, state of health can be correlated with impedance data. The "dead cell", which 
presumably failed because of decreased active zinc content and decreased active zinc surface area, 
had a significantly higher low frequency impedance. 

3. Semi-empirical equivalent circuit models were able to fit the experimental data. It is 
recommended that models incorporating the fundamental coupled differential equations describing 
the electrodes be solved in future work. Such an approach will be more difficult, and require 
significant time and resouces to initiate, but the results will be more fundamentally satisfying. 
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Table 1. Equivalent Circuit Parameters. 



Silver 

Nickel 

Cadmium 

aC 

212 F 

0.021 F 

13.25 F 

R 

f 

0.083 Q 

0.0062 n 

0.044 Cl 

A 

2.2x10^ Q s “ 

2.4xl0 5 ns m 

1.1x10 s Cl s 

m 

-1.1 

-1.7 

-1.9 
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